How complex is the quantum motion? 
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In classical mechanics the complexity of a dynamical system is characterized by the rate of local 
exponential instability which effaces the memory of initial conditions and leads to practical irre- 
versibility. In striking contrast, quantum mechanics appears to exhibit strong memory of the initial 
state. Here we introduce a notion of complexity for a quantum system and relate it to its stability 
and reversibility properties. 

PACS numbers: 05.45.Mt, 03.65.Sq, 05.45.Pq 



The question of how complex is quantum motion is 
of fundamental importance with deep connections to en- 
tanglement and decoherence. However, our knowledge 
of the relations between complexity, dynamical stabil- 
ity, reversibility and chaos is far from being satisfactory, 
sometimes even confusing and a clearer understanding is 
necessary. 

To this end let us first consider classical motion where 
things arc quite well settled. Classical complex systems 
are characterized by positive Lyapunov exponent i.e. by 
local exponential instability. They have positive algo- 
rithmic complexity and, in terms of the symbolic dy- 
namical description, almost all orbits are random and 
unpredictable 

In spite of many efforts Q , the problem of character- 
izing the complexity of a quantum system is still open. 
Indeed the above notion of complexity cannot be trans- 
ferred, sic et simpliciter, to quantum mechanics, where 
there is no notion of trajectories. Still, a comparison 
between classical and quantum dynamics can be made 
by studying the evolution in time of the classical and 
quantum phase space distributions, both ruled by linear 
equations. 

First investigations focused on reversibility, namely on 
the propagation of round-off errors in numerical simula- 
tions [3| . Strong and impressive evidence has been gath- 
ered that the quantum evolution is very stable, in sharp 
contrast with classical dynamics in which the extreme 
sensitivity to initial conditions, which is the very essence 
of classical chaos, leads to a rapid loss of memory. 

Later on, a different approach focussed on the stabil- 
ity properties of motion under small variations of system 
parameters. This approach does not rise any difficulty 
in the classical context since exponentially unstable sys- 
tems exhibit the same rate of exponential instability by 
slightly changing the initial conditions with fixed param- 
eters or by changing parameters with fixed initial con- 
ditions. On the other hand the advantage of the latter 
approach is that it can be applied to phase space distri- 
butions. Here one computes the so-called fidelity 0, Q , 
defined as the overlap between two distributions evolving 



under two slightly different Hamiltonians. It is tempt- 
ing to connect the behavior of fidelity to the regular or 
chaotic behavior of quantum motion. Indeed the orig- 
inal expectation, which seemed quite natural, was that 
fidelity should remain close to one at all times for inte- 
grable systems and fall down exponentially for chaotic 
systems However, this expectation is not fulfilled. 

For the purpose of the present paper it is necessary 
here to make clear the following. The motivation for the 
introduction of fidelity, namely the suggestion to analyze 
the stability of motion by perturbing the Hamiltonian 
rather than the initial state, originated from the obser- 
vation that in quantum mechanics, due to the unitary 
evolution, the scalar product (i[>(t)\i/j' {t)} of two initially 
close states |^>(0)) and |^'(0)) does not change in time. 
However two points must be stressed: i) the classical 
evolution of phase space density is also unitary and lin- 
ear and therefore the overlap of two initially close phase 
space distributions does not change with time in classical 
mechanics as well; ii) for classically chaotic quantum sys- 
tems, the fidelity decay, depending on the perturbation 
strength, can be Gaussian or exponential. A power-law 
decay is also possible in the quantum diffusive regime [1] . 
Furthermore, for integrable systems the fidelity decay can 
be faster than for chaotic systems [j| . On the other hand, 
even in classical mechanics the fidelity decay does not 
clearly distinguish between chaotic and integrable sys- 
tems. In short, fidelity is not a good quantity to charac- 
terize the complexity of motion, neither in quantum nor 
in classical mechanics. 

In this paper, we propose the number of harmonics of 
the Wigner function as a suitable measure of the com- 
plexity of a quantum state. We recall that in classical me- 
chanics the number of harmonics of the classical distribu- 
tion function in phase space grows linearly for integrable 
systems and exponentially for chaotic systems, with the 
growth rate related to the rate of local exponential insta- 
bility of classical motion Q . Thus the growth rate of the 
number of harmonics is a measure of classical complex- 
ity. Since the phase space approach can be equally used 
for both classical and quantum mechanics, the number 
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of harmonics of the Wigner function appears as the cor- 
rect quantity to measure the complexity of a quantum 
state. In what follows, we examine the behaviour of this 
quantity and its relation to fidelity and reversibility prop- 
erties. Moreover, we show that the number of harmonics 
can be used to detect the transition from integrability 
to quantum chaos. A detailed derivation of the results 
discussed in this paper can be found in Ref. [Io| . 

The Wigner Function. - Let us consider a generic 
nonlinear system which exhibits a transition from quasi- 
integrable to chaotic behavior as the strength of the 
nonlincarity is increased. More precisely we consider 
the Hamiltonian operator H = H(at,a;t) = H^'in — 
t) with a time-independent unperturbed 
part H^ ' with a discrete energy spectrum bounded from 
below. Here d\a are the bosonic, creation-annihilation 
operators and [a, at] = 1 . 

We will use the method of c-number a-phase space bor- 
rowed from quantum optics (see for example This 
method is basically built upon the basis of the coher- 
ent states | a) which are defined by the eigenvalue prob- 
lem a\a) — -^\ a )i where a is a complex variable inde- 
pendent of the effective Planck's constant h. An arbi- 
trary coherent state is obtained from the ground state, 
\a) = D |0), with the help of the unitary displace- 

ment operator D (A) = cxp(Afit — A*d). The Wigner 
function W in the a-phase plane is related to the density 
operator p as follows: 



W(a*,a;t) 
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(1) 

where the integration runs over the complex 77-plane. 

The harmonic's amplitudes W m of the Wigner function 
are given by the expansion 
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W(a*,a;t) = - V W m (I;t), 



(2) 



where a = \/l e~ lB , with (1,9) action-angle variables. 

The fidelity. - Following the approach developed in 
Rcf. I12L we consider now the forward evolution 



p(t) = U{t)p(Q)U\t) 



(3) 



of an initial (generally mixed) state /5(0) up to some time 
t = T. A perturbation P(£) is then applied at this time, 
with perturbation strength £. For our purposes, it will 
be sufficient to consider unitary perturbations P(£) = 
e~ l ^ v , where V is a Hermitian operator. The perturbed 
state 



~p(t,o = p(Op(t)PHO 



(4) 



is then evolved backward, with the same Hamiltonian, 
for the time T, thus obtaining the final state 



fi(0\T,O = U*(T)fi(T,Z)U(T). 



Finally, we consider the distance between the reversed 
p(0\T, £) and the initial p(0) state, as measured by the 
Peres fidelity [H 

m T) = ■fr[p(Q|r,flp(Q)] = ^feOOT (6) 



Tr[p2(0)] 



Tr[p 2 (T)] 



This quantity is bounded in the interval [0, 1] and the 
distance between the initial and the time-reversed state 
is small when F(£;T) is close to one. In particular, 
F(£; T) = 1 when the two states coincide. The last equal- 
ity in ([6]) is a consequence of the unitary time evolution. 

The Peres fidelity ^ can be expressed in terms of the 
Wigner function as 



F&T) 



J d 2 a W (a*, a; 0) W (a*, a; 0\T, f) 

J d 2 aW 2 (a*,a;0) 
/ d 2 a W (a*, a; T) W (a*, a; T, £) 
/ d 2 aW 2 (a*,a;T) ' 



(7) 



The advantage of this representation is that it remains 
valid in the classical case when the Wigner function re- 
duces to the classical phase-space distribution function 
W c (a*,a;t). 

Growth of the number of harmonics.- As is well known, 
the paramount property of the classical dynamical chaos 
is the exponentially fast structuring of the system's phase 
space on finer and finer scales. In particular, the num- 
ber M.(i) of harmonics (Fourier components) that sig- 
nificantly contribute in the expansion ([2]) of the classical 
phase space distribution W C7 grows exponentially in time 
for chaotic motion. The crucial point is that only in quan- 
tum mechanics the number of harmonics of the Wigner 
function is directly related to the expectation value of 
physical observables. Therefore in quantum mechanics 
an exponential growth of number of harmonics is not al- 
lowed in general 0, [3, 14 1 . 

For an explicit numerical evaluation of the number of 
harmonics we will focus on the quantity 



(m 2 )t 
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The quantity \J (m 2 ) t provides an estimate of the number 
M. (t) of harmonics developed by the time t and therefore 
of the complexity of the Wigner function at time t. 

Relation between fidelity and number of harmonics. - 
This relation takes a very simple form when the pertur- 
bation at the reversal time t = T is a rotation of the 
quantum phase-space distribution W(a* ,a;t). Consider 
the unitary transformation (|3|) with the perturbation op- 
erator V = n and the rotation angle £. In this case we 
obtain [l(| 



F(fct) = l-2 



E+oo 
m=- 
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(£m/2) f Q °° dI\W m (I;t)\* 



(5) 



El=-o Io°°dI\W m (I;t^ 



(9) 



3 



The lowest-order ^-expansion of this equation reads 



1 



F(C;t)«l-^ 
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(10) 



Notice that relations ([9]) and (fTO)) between distance 
-F and number of harmonics can be applied to classical 
dynamics, provided that the harmonics of the classical 
distribution function W c instead of those of the Wigncr 
function W are used. 

Illustrative example-. Let us consider the kicked 
tic oscillator model, defined by the Hamiltonian 
IE El 

H = hujQn + h 2 h 2 - \fhg(t)(a + (11) 
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where g(t) = go^2 s S(t — s), h = a^a, [a,aJ] = 1. In 
our units, the time and parameters h, loq as well as the 
strength go of the driving force arc dimcnsionlcss. The 
period of the driving force is set to one. The correspond- 
ing classical Hamiltonian function, expressed in terms of 
complex canonical variables a,ia*, is given by 



H c = uj \a\ 2 + H 4 - g(t){a* + a) . 



(12) 



We choose the initial state to be an isotropic mixture 
of coherent states: p(0) = j d?aV(\a\ 2 )\a)(a\, where 
V(l = \a\ 2 ) = e _/ / A . Correspondingly, the initial 
Wigncr function is isotropic and Gaussian, W(a* , a; 0) oc 



e ~ a+h/2 _ Only the m = harmonic is excited, so that 
(m 2 )t=o = 0. The particular case A = corresponds 
to the pure ground state, which occupies the minimal 
quantum cell with the area h/2 . 

The classical dynamics of model (fTTj) becomes chaotic 
(with negligible stability islands) when the perturbation 
strength go > 1. In this regime the mean action grows 
diffusively with the diffusion coefficient D w g 2 . 

Chaotic regime. - Wc now compare, in the chaotic 
regime, the evolution in time of (m 2 ) t for quantum 
and classical dynamics. For this purpose, we solve 
both the quantum and the classical Liouville equation. 

In the latter case, the initial phase space distribution 

_id 2 

W c (a*,a;0) oc e — s~~ has size S which coincides, for a 
given value of H, with the size h/2 of the Wigner func- 
tion corresponding to the initial quantum ground state 
p(0) = |0)(0|. The quantum to classical transition is 
then explored by keeping S constant and considering, for 
smaller and smaller values of h, initial incoherent mix- 
tures of size S = A + h/2. The results are shown in Fig. [I] 
The exponential increase of (m 2 } t takes place only up to 
the Ehrenfest time scale oc In h 11511 , consistently with 
the findings reported in Refs. [^, |T3|. After that time, 
a much slower power-law increase follows. Namely the 
number of harmonics .M(i) ~ \J (m 2 ) increases linearly 
for the pure state case, where (m 2 )t ~ oc t 2 due to 
diffusive growth of the mean action, and slower than lin- 
early for mixtures [ltj| . This gro wth eventually saturates 
due to quantum localization [13j of diffusive motion. 




20 25 



FIG. 1: (color online) Root-mean-square radius (m }t of the 
distribution of harmonics versus time t, at wo = 1, ga = 
1.5, S = 0.5. Squares, diamonds and triangles correspond to 
h = 0.01,0.1 and 1. In this latter case, p(0) = |0)(0|. Empty 
circles refer to classical dynamics and the dashed line fits these 
data. 



From Eq. (|10p we can estimate a critical perturba- 
tion strength £ C (T) w ^/2/(m 2 )r, such that the fidelity 
F(£;T) remains close to 1 after the backward evolution 
as long as £ -C £c(T), whereas reversibility is lost when 
£ ^ £c(T). This statement is illustrated in Fig. [5] for the 
number of harmonics, estimated by */ (m 2 ) t . Therefore, 
we establish a direct connection between complexity of 
phase space distribution and degree of reversibility of mo- 
tion. Due to the strikingly different growth in time of the 
number of harmonics for classical and quantum chaotic 
motion, £ C (T) drops exponentially with T in the classical 
case and at most linearly in the quantum case (after the 
Ehrenfest time scale). Therefore, our analysis explains 
the numerically observed Q much weaker sensitivity of 
quantum dynamics to perturbations as compared to clas- 
sical dynamics. 

Note that, due to the exponential proliferation of the 
number of harmonics in classical mechanics, formula (|10p 
is only valid up to a time logarithmically short in the 
perturbation strength £. After that time, due to diffu- 
sive growth of the mean action, the decay of the fidelity 
F turns from exponential to power law [l8(, F oc l/(Dt), 
while the number of harmonics still grows exponen- 
tially with time and correctly describes the complexity 
of chaotic motion. 

Crossover from integrability to quantum chaos. - It 
is known that in the integrable regime the number of 
harmonics, computed in the action-angle representation, 
grows linearly with time Q and that nearby orbits sepa- 
rate linearly fast [lj)[ . We show (see the inset of Fig. [3]) 
that the linear growth of J (m 2 ) also takes place in 
quantum mechanics, up to the Heisenberg time scale 
tn oc h~ x . The strikingly different behavior of the num- 
ber of harmonics of the Wigner function in the integrable 
and chaotic regime suggests that this quantity can be 
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FIG. 2: (color online) Reversibility properties of quantum 
dynamics. The backward evolution starts at the reversal 
time T = 50. We show \J (m 2 )t, for different values of the 
perturbation parameter: from bottom to top, £ = £ C (T) x 
cxp (—1/2), I = 8, . . . , 1, / = (thick black curve marked by 
the closed circle), and I = —1, . . . , —6, at u>o = 1, h = 1, 
9o = 
each 



ncntial instability. Wc would like to stress that in relation 
to other, very interesting, proposed measures of quantum 
complexity, such as quantum dynamical entropies [H , our 
quantity is very convenient for numerical investigations. 
It becomes therefore possible to investigate complexity as 
a function of the effective Planck's constant. Moreover, 
the above outlined phase-space approach is quite general 
and can be readily extended to systems with arbitrary 
number of degrees of freedom, including qubit systems, 
whose Hamiltonian can be expressed in terms of a set 
of bosonic creation-annihilation operators. Therefore, in 
many-body systems the (number of) harmonics of the 
Wigncr function could shed some light on the connection 
between complexity and entanglement, a fundamental is- 
sue of great relevance for the prospects of quantum in- 
formation science. Moreover, we believe that the depen- 
dence of the number of harmonics on control parameters 
could be used not only to investigate the integrability 
to quantum chaos crossover but also to detect quantum 
phase transitions. 
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FIG. 3: (color online) Dependence of {m 2 )t on the pertur- 
bation strength go at time t — 3, for h = 0.01, loo = 1, 
A = 0. Inset: time evolution of %J (m 2 )t in the integrable 
regime, at go — 1, u>o = 1, 8 = 0.5 and, from bottom to top, 
h=l, 0.1, 0.05, 0.02, 0.01, 0.005. 

used to detect the transition to quantum chaos. We 
show in Fig. [3] that (to 2 ), computed at a given time t, 
exhibits a sharp increase when the perturbation param- 
eter <?o > 0.5. From this figure we can conclude that the 
crossover from integrability to chaos takes place in the 
region 0.5 < go < 0.7. 

Discussion. - To summarize, we have shown that the 
number of harmonics of the Wigner function is a suitable 
measure of complexity of a quantum state, in that this 
quantity is directly related to the reversibility properties 
of quantum motion and, at the classical limit, reproduces 
the well-known notion of complexity based on local expo- 
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